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Abstract 

Let G = (Cti,Ct2) be a Gaussian vector in R 2 with EG1G2 7^ 0. 
Let ci,C2 £ R . A necessary and sufficient condition for G = ((Gi + 
c\a) 2 , {G2 + C2«) 2 ) to be infinitely divisible for all a G R 1 is that 

r M >-r^>o vi<i/j<2. (0.1) 

In this paper we show that when (jO.ip does not hold there exists an 
< ceo < 00 such that G = ((G\ + C\a) 2 , (G2 + c^a) 2 ) is infinitely 
divisible for all |q| < ao but not for any \a\ > oq. 

1 Introduction 

Let 7] = (rjx, ... , rj n ) be an R n valued Gaussian random variable, rj is said to 
have infinitely divisible squares if rf := (rj 2 , . . . ,t] 2 ) is infinitely divisible, i.e. 
for any r we can find an R n valued random vector Z r such that 

rj 2 l = Z r>1 H h Z r>r , (1.1) 

where {Z r j}, j = 1, . . . ,r are independent identically distributed copies of 
Z r . We express this by saying that rf is infinitely divisible. 
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Paul Levy proposed the problem of characterizing which Gaussian vectors 
have infinitely divisible squares. It is easy to see that a single Gaussian 
random variable has infinitely divisible squares. However, even for vectors 
in R 2 this is a difficult problem. It seems that Levy incorrectly conjectured 
that not all Gaussian vectors in R 2 have infinitely divisible squares. If he 
had said R 3 his conjecture would have been correct. 

Levy's problem was solved by Griffiths and Bapapt [Tj, [7J, (see also [HI 
Theorem 13.2.1]). 

Theorem 1.1 Let G = (G\, . . . , G n ) be a mean zero Gaussian random vari- 
able with strictly positive definite covariance matrixT = {Tjj} = {E(GiGj)}. 
Then G 2 is infinitely divisible if and only if there exists a signature matrix 
Af such that 

A/T'W is an M matrix. (1.2) 

We need to define the different types of matrices that appear in this 
theorem. Let A = {a-i,j}i<i,j< n be an n x n matrix. We call A a positive 
matrix and write A > if djj > for all i,j. We say that A has positive 
row sums if J2]=i a i,j > for all 1 < i < n. 

The matrix A is said to be an M matrix if 

(1) a it j < for all i ^ j. 

(2) A is nonsingular and A^ 1 > 0. 

A matrix is called a signature matrix if its off-diagonal entries are all 
zero and its diagonal entries are either one or minus one. The role of the 
signature matrix is easy to understand. It simply accounts for the fact that 
if G has an infinitely divisible square, then so does (eiGi, . . . , e n G n ) for any 
choice of q = ±1, i = 1, . . . , n. Therefore, if ( 11.21) holds for M with diagonal 
elements ni, . . . , n n 

(MY- 1 Af)' 1 = ATT AT > (1.3) 

since the inverse of an M matrix is positive. Thus (n\Gi, . . . ,n n G n ) has a 
positive covariance matrix and its inverse is an M matrix. (For this reason, 
in studying mean zero Gaussian vectors with infinitely divisible squares one 
can restrict ones attention to vectors with positive covariance.) 

The natural next step was to characterize Gaussian processes with in- 
finitely divisible squares which do not have mean zero. We set rji = G{ + q, 
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EGi = 0, % = 1, , . . . , n. Let T be the covariance matrix of (G\, . . . , G n ) and 
set 

c := (ci, . . . ,c n ). (1.4) 

Set 

G + c:= (Gi + ci,...,G r n + c n ) (1.5) 

and 

(G + c) 2 := ((d + Cl ) 2 , . . . , (G n + c n f). (1.6) 

Results about the infinite divisibility of (G + c) 2 when c\ = ■ ■ ■ = c n , 
are given in the work of N. Eisenbaum [21 E] and then in joint work by 
Eisenbaum and H. Kaspi [I], as a by product of their characterization of 
Gaussian processes with a covariance that is the 0-potential density of a 
symmetric Markov process. We point out later in this Introduction how 
Gaussian vectors with infinitely divisible squares are related to the local times 
of the Markov chain that is determined by the covariance of the Gaussian 
vector. It is this connection between Gaussian vectors with infinitely divisible 
squares and the local times of Markov chains, and more generally, between 
Gaussian processes with infinitely divisible squares and the local times of 
Markov processes, that enhances our interest in the question of characterizing 
Gaussian vectors with infinitely divisible squares. 

Some of the results in [2J [3], H] are presented and expanded in 0, Chapter 
13]. The following theorem is taken from [9], Theorem 13.3.1 and Lemma 
13.3.2]. 

Theorem 1.2 Let G = (G\, . . . , G n ) be a mean zero Gaussian random vari- 
able with strictly positive definite covariance matrixY = {T it j} = {E(GiGj)}. 
Let 1 = (1, . . . , 1) G R n . Then the following are equivalent: 

(1) (G + la) has infinitely divisible squares for all a e R 1 ; 

(2) For £ = N(0,b 2 ) independent of G , (G\ + £,..., G n + £, £) has infinitely 
divisible squares for some b ^ 0. Furthermore, if this holds for some 
6/0, it holds for all b 6 R 1 , with N(0, 0) = 0. 

(3) T^ 1 is an M matrix with positive row sums. 

In [8], Theorem 11.21 is generalized so that the mean of the components of 
G + c in (jl.5p need not be the same. In this generalization certain trivial 
cases spoil the simplicity of the final result. We avoid them by requiring that 
the covariance matrix of the Gaussian process is irreducible. 
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Theorem 1.3 Let G = (G\, . . . , G n ) be a mean zero Gaussian random vari- 
able with irreducible strictly positive definite covariance matrix V = {Tij} = 
{E{GiGj)}. Let c = (ci, . . . , c n ) G R n , and let C be a diagonal matrix 
with Ci = Cij, 1 < i < n . Then the following are equivalent: 

(1) G + ca has infinitely divisible squares for all a G R 1 ; 

(2) For £ = iV(0,& 2 ) independent of G, (G x + c x £, ...,G n + c n £,£) has 
infinitely divisible squares for some 6^0. Furthermore, if this holds 
for some b ^ 0, it holds for all b G R 1 ; 

(3) CT~ l C is an M matrix with positive row sums. 

By definition, when (G + c) 2 is infinitely divisible, it can be written as 
in fll.ip as a sum of r independent identically distributed random variables, 
for all r > 1. Based on the work of Eisenbaum and Kaspi mentioned above 
and the joint paper [5] we can actually describe the decomposition. We give 
a rough description here. For details see [21 [HH] and [91 Chapter 13]. 

Assume that (1), (2) and (3) of Theorem 11.31 hold. Let 



Let T c denote the covariance matrix of G/c. Theorem 11.21 holds for G/c 
and r c , so T^ 1 is an M matrix with positive row sums. To be specific let 
G/c G R n . Set S = {1, . . . ,n}. By [9, Theorem 13.1.2], T c is the 0-potential 
density of a strongly symmetric transient Borel right process, say X, on S. 
We show in the proof of [9l Theorem 13.3.1] that we can find a strongly 
symmetric recurrent Borel right process Y on S U {0} with P X (T < oo) > 
for all x G S such that X is the process obtained by killing Y the first time 
it hits 0. Let L x t = {L x ; t G R+, x G S U {0}} denote the local time of Y. It 
follows from the generalized second Ray-Knight Theorem in [5], see also [9j 
Theorem 8.2.2] that under P° x P G , 



for all t G R + , where r(t) = inf{s > 0|L° > t}, the inverse local time at zero, 
and Y and G are independent. Consequently 




(1.7) 




(1.8) 




(1.9) 
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for all a G R 1 . (We can extend a from R + to R 1 because G is symmetric.) 
{ c x^T(a 2 /2)i x e ^} anc ^ ^ £ 5 1 } are independent. G 2 is infinitely 

divisible and for all integers r > 1 

c2 L T ( a 2/ 2 ) = C 2 I/ r ( a 2/( 2r )) j i + 1" C 2 -^ T ( Q ,2/( 2r )), r (1-10) 

where {-^ T ( a 2 / ( 2r , ) ) J }, j = 1, . . . ,r are independent. 

Note that in (11. 9p we identify the components of the decomposition of 
{(G x + c x a) 2 ; x E S} that mark it as infinitely divisible. 

In Theorem 11.31 we have necessary and sufficient conditions for ((Gi + 
Cia) 2 , (Cr 2 + c 2 a) 2 ) to be infinitely divisible for all a 6 fi 1 . There remains 
the question, can ((Gi + c\a) 2 , (G 2 + c 2 a) 2 ) have infinitely divisible squares 
for some a > but not for all a G -R 1 ? When we began to investigate 
this question we hoped that such points a do not exist. This would have 
finished off the problem of characterizing Gaussian random variables with 
infinitely divisible squares and, more significantly, by (11.91) . would show that 
when a Gaussian random variable with non-zero mean has infinitely divisible 
squares, it decomposes into the sum of two independent random variables. 
The Gaussian random variable itself minus its mean, and the local time of a 
related Markov process. This would be a very neat result indeed, but it is 
not true. 

For all Gaussian random variables (G\, G 2 ) in R 2 and all c\, c 2 G -R 1 define 

Q 2 ( Cl , c 2 , a) := ((d + Cl a) 2 , (G 2 + c 2 a) 2 ). (1.11) 

It follows from Theorem ll.3[ (for details see [8j Corollary 1.3, 4.]), that 
G 2 {ci, c 2 , a) has infinitely divisible squares for all a G -R 1 if and only if 

ri,i>-r lj2 and r 2j2 >^r li2 . (1.12) 

C 2 Ci 

If (11.121) does not hold, we call < «o < oo a critical point for the infinite 
divisibility of G 2 (ci, c 2 , a) if G 2 (ci, c 2 , a) is infinitely divisible for all \a\ < a , 
and is not infinitely divisible for any \a\ > a^. In this paper we prove the 
following theorem: 

Theorem 1.4 For all Gaussian random variables (Gi,G 2 ) in R 2 and all 
(c%, c 2 ) G R 2 for which U.lty) does not hold, G 2 (ci, c 2 , a) has a critical point. 
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Note that in Theorem 11.41 we consider all (ci,c 2 ) G R 2 . It follows from 
(j!.12p that when EG1G2 > 0, then Q 2 (ci, c 2 , a) has infinitely divisible squares 
for all a G R l only if cic 2 > 0. Nevertheless, by Theorem 11.41 even when 
C1C2 < 0, (Gi + cia,G2 + c 2 a) does have infinitely divisible squares for |a| 
sufficiently small. 

To conclude this Introduction we explain how we approach the problem 
of showing that [G + ca) 2 is infinitely divisible for all a G R 1 or only for 
some a G R 1 . Since we can only prove Theorem 11.41 for Gaussian random 
variables in R 2 we stick to this case, although a similar analysis applies to 
R n valued Gaussian random variables. 

Let T be the covariance matrix of G and 

f := (I + TA)- 1 T = (T- 1 + A)~\ (1.13) 

where 

(Li4) 

Consider the Laplace transform of ((Gi + Ci«) 2 , (G 2 + c 2 a) 2 ), 

3 -(Ai(Gi+cia) 2 +A2(G2+C2a) 2 )/2 N 



1 I o? ( 2 

exp I — — I c 2 Ai + C2A2 - CiCjXiXjT^ 



(det(/ + rA))V2 ^1 2 



(See [3 Lemma 5.2.1].) Set Ai = t(l — si) and A 2 = t(l — s 2 ), < s%, s 2 < 1 
and write (jl.150 as 



exp 



([/ (si, s 2 , t, r) + « 2 v^ 8 2 , t, ci, c 2 , r)) (1.16) 



where C/ := f/(s 1; s 2 , t, T) = -1/2 log(det (J + TA)). Note that 
exp {U (si, s 2 , t, T)) is the Laplace transform of (G 2 , G 2 ), with the change of 
variables Ai = t(l — si) and A 2 = t(l — s 2 ). It is easy to see from Theorem ll.il 
that all two dimensional Gaussian random variables are infinitely divisible. 
Therefore, for all t sufficiently large, all the coefficients of the power series 
expansion of U in Si and s 2 are positive, except for the constant term. This is 
a necessary and sufficient condition for a function to be the Laplace transform 
of an infinitely divisible random variable. See e.g. [91 Lemma 13.2.2]. 

Now, suppose that for all t sufficiently large, V (si, s 2 , t, ci, c 2 , f ) has all 
the coefficients of its power series expansion in si and s 2 positive, except 
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for the constant term. Then the right-hand side of (I1.15P is the Laplace 
transform of two independent infinitely divisible random variables. It is 
completely obvious that this holds for all a G R 1 . 

On the other hand suppose that for all t sufficiently large, the power series 
expansion of V (sx, s 2 , t, ci, c 2 , V) in si and s 2 has even one negative coeffi- 
cient, besides the coefficient of the constant term. Then for all a sufficiently 
large, (11.161) is not the Laplace transform of an infinitely divisible random 
variable. In other words ((Gi+cia) 2 , (G 2 +c 2 a) 2 ) is not infinitely divisible for 
all a G R 1 . But it may be infinitely divisible if a is small, since the positive 
coefficients of U may be greater than or equal to a 2 times the corresponding 
negative coefficients of V. Clearly, if this is true for some |a| = «o > 0, then 
it is true for all \a\ < a . 

The preceding paragraph explains how we prove Theorem 11.41 We con- 
sider vectors ((Gi + cia) 2 , {G2 + 02a) 2 ) that are not infinitely divisible for all 
a G R 1 , (this is easy to do using (I1.12p ). and show that for |a| sufficiently 
small the coefficients in the power series expansion of 

U (si, s 2 , t, T) + a 2 V (si, s 2 , t, ci, c 2 , f) (1.17) 

in Si and s 2 are positive, except for the constant term. Our proof only uses 
elementary mathematics, although it is quite long and complicated. In the 
course of the proof we show that the coefficients of the power series expansion 
of U in si and s 2 are positive, except for the constant term. This provides 
a direct elementary proof of the fact that the Gaussian random variable 
(Gi,G 2 ) always has infinitely divisible squares. 

As we have just stated, and as the reader will see, the proof of Theorem 
II. 41 is long and complicated. So far we have not been able to extend it to apply 
to Gaussian random variables in R 3 . One hopes for a more sophisticated and 
much shorter proof of Theorem 11.41 that doesn't depend on the dimension of 
the Gaussian random variable. 

2 Gaussian squares in R 2 and their Laplace 
transforms 

Let G = (Gi, G2) be a mean zero Gaussian process with covariance matrix 

r=(;i) P.i) 
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where ab = d + 1 > 1, and let G + c := (G\ + c\, G 2 + c 2 ), c\, c 2 G i? 1 . Note 
that 

detr = rf, (2.2) 



and 



Let 



Then 



and 



-■4 v: <»> 



A = ' A I ' <-) 
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r- 1 + A = i f 6 + f Al -1 ) (2.5) 



r := (/ + rA) _1 r = (r _1 + A) -1 (2.6) 

' b + d\ 1 -1 y 1 

— 1 a + d\ 2 J 
1 f X 2 d + a 1 \ 



H(a, b, X 1 ,X 2 ) \ 1 X±d + b 



where 



F(a, 6, Ai, A 2 ) = 1 + aAi + 6A 2 + dX t X 2 = d det (r -1 + A) . (2.7) 

(We use repeatedly the fact that ab = d + 1.) Note that by (12. 2p we have 
that 



det (/ + TA) = det (T (T^ + AJJ = ddet (r^ 1 + AJ = H(a, b, A 1; A 2 ). 

(2.8) 

Lemma 2.1 

E G (e" (Al(Gl+Cl)2+A2(G2+C2)2)/2 ) (2.9) 
I exp ^ Ai + c 2 A 2 + (ci& + c 2 a - 2c x c 2 ) AiA 2 ^ 



(#(a,Mi,A 2 ))V2 *y 2H(a,b,X 1 ,X 2 ) 
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Proof By P Lemma 5.2.1] 

E G ^ e -(Ai(Gi+ C i) 2 +A2(G 2 + C2 ) 2 )/2^ ^2.10) 




■exp I --( c\\x + c 2 A 2 



(H(a, b, Ai, A 2 )) 1/2 

cjXj (\ 2 d + a)+ 2 Cl c 2 AiA 2 + c|A|(Aij + 6) \ \ 
tf(a,Mi,A 2 ) JJ ' 

A simple computation shows that 

(c?A 1 + 4\ 2 )#(a,&,A 1 ,A 2 ) (2.11) 
- (c 2 A 2 (A 2 d + a) + 2cic 2 AiA 2 + c^Aid + &)) 
= c\\i + c\\ 2 + {c\b + c\a - 2cic 2 ) AiA 2 , 

from which we get (12.91) . □ 

The term l/(H(a, b, X 1 , A 2 )) 1//2 is the Laplace transform of [G\^G 2 ^)j2 
and by Corollary 1.1, 2.] it is the Laplace transform of an infinitely 
divisible random variable. The exponential term may or may not be a Laplace 
transform. In fact, by (11.121) . we know it is the Laplace transform of an 
infinitely divisible random variable, for all {cia,c 2 a}, for all a G R 1 , if and 
only if 

a> — >0 and b > — > 0. (2.12) 

c 2 Ci 

To prove Theorem 11.41 we must show that when (12.121) does not hold, there 
exists an < «o < oo such that (12. 9p is the Laplace transform of an infinitely 
divisible random variable when c\ and c 2 are replaced by C\a and c 2 a for any 
| a | < ao- Actually, as we see in Section [HI the general result follows from the 
consideration of three cases, 

1. Ci = c 2 := (c,c); 

2. ci = -c 2 := (c, -c) 

3. (c,0). 

This is because if ci 7^ c 2 and neither of them is zero, we can replace 
(Gi,G' 2 ) by (Gi/|ci|, G 2 /|c 2 |). Clearly, in this case, if Theorem 11.41 holds 
for (Gi/|ci|,G 2 /|c 2 |) it holds for (G 1 ,G 2 ). 

In these three cases the numerator of the fraction in the exponential term 
on the right-hand side of (12.91) is 
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1. c 2 ((a + 6- 2)A 1 A 2 + A 1 + A 2 ); 

2. c 2 ((a + 6 + 2)A 1 A 2 + A 1 + A 2 ) 

3. c 2 (6AiA 2 + Ai). 
Set 

1 = a + b-2 and p = a + b + 2. (2.13) 

Note unless detT = 0, ab > 1. Since Theorem 11.41 obviously holds when 
det T = 0, we can exclude this case from further consideration. Thus we 
always have 7 > 0. 



3 Power series expansion of the logarithm of 
the Laplace transform of ((Gi + c) 2 , (G2 + c) 2 
when EG1G2 = 1 

Bapat's proof of Theorem 11.11 involves the analysis of a certain power series 
expansion of the logarithm of the Laplace transform. We need a similar, but 
more delicate, analysis. (See Lemma [4.11 below) . 

Using (I2.9p and the remarks following Lemma 12.11 we can write 

E G (V (Al(Gl+c)2+A2(G2+c)2)/2 ) (3.1) 

/ 1 1 u( u \ \ \ \ ( c2 (7AiA 2 + Ax + A 2 )^ 
= exp --logiJ(a, 0, Ai,A 2 ) exp 1 



2 ^ ^ ' " V 2H(a,b,X h X 2 ) 

exp ( -(P(a, b, Ai, A 2 ) + c 2 Q(a, b, A x , A 2 ) 



.2 

Since ab = d + 1, (recall that d > 0), we have 

a + b-{d + 2) = a+ ^ + ^ - (d + 2) (3.2) 

a 

= - (a 2 - (d + 2)a + (d+ 1] 

(X 

= I( a -(d + l))(a-l). 
a 

Thus a + b — (d + 2) < if and only if 1 < a < d + 1, which in view of 
ab = d + 1 is equivalent to 1 < b < d + 1. Consequently (I2.12p holds if and 
only if a + b - (d + 2) < 0. 
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Therefore, to show that ((Gj + c) 2 , (G 2 + c) 2 ) is infinitely divisible, for 
some, but not for all, c > 0, we must consider a, b > such that 

( ■= a + b- (d + 2) > 0. (3.3) 

In the rest of this paper we assume that (13.31) holds. 

Let Ai = t(l — si) and A2 = t(l — s 2 ), < si, s 2 < 1. We consider P and 
Q as functions of and s\, s 2 , t, and write 

P(s 1 ,s 2 ,t) :=P(a, b, Ai,A 2 ), g(si,s 2 ,t) :=Q(a, 6, Ai,A 2 ). (3.4) 

We expand these in a power series in si, s 2 . 

00 00 

P(si,s 2 ,t) = Pi,k(t)44, Q(si,s 2 ,t) = £ Qi,fc(*M4 (3-5) 
j,fc=o i,fe=o 

and set 

J?(si, s 2 , t, c) = P(si, s 2 , t) + c 2 Q(si, s 2 , t, c). (3.6) 

Consequently 

00 

fl(si,s 2 ,t,c)= Rj,k(^ c )44 ( 3 - 7 ) 

j,fc=0 

with 

R j>k {t,c) = P jik {t)+c 2 Q j!k {t). (3.8) 
In this section we obtain explicit expressions for Pj jk (t), Qj,k(t). We write 

H(a,b,X 1 ,X 2 ) = 1 + aAi + 6A 2 + rfAiA 2 (3.9) 
= 1 + at + bt + dt 2 - (at + dt 2 )si - (bt + dt 2 )s 2 + dt 2 SiS 2 
= (1 + at + bt + dt 2 )(l — asi — [3s 2 + psis 2 ) 
= (l + at + bt + dt 2 )(l - asi - (3s 2 + 0al3siS 2 ) 

where 

at + dt 2 n bt + dt 2 dt 2 

a = ^, — T—r^, ,.9 > P = ; , : , u , rz? > P 



1 + at + bt + dt 2 1 + at + bt + dt 2 1 + at + bt + dt 2 

(3.10) 

and 

(l + a 1 - 1 + at + bt + dt 2 \ p , 011 . 
^ = P 1 + at + bt + dt 2 '= 6- (3 - U) 
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Note that 

d- 1 1 
1 ~ 6 = 1 + d- 1 + at + bt + dt 2 ~ ~dPP' (3 ' 12) 
Using these definitions we have 

P(a, b, Ai,A 2 ) = - log if (a, 6, Ai,A 2 ) (3.13) 
= - log(l + at + bt + rft 2 ) - log(l - asi - (3s 2 + 9af3s 1 s 2 ) 



and 



Q(a,6, Ai, A 2 ) (3.14) 
7A1A2 + Ai + A 2 
H(a, b, Ai, A 2 ) 

1 7 t 2 (l-s 1 )(l-s 2 )+t(2-s 1 -s 2 ) 



(1 + at + bt + <it 2 ) 1 — asi — (3s 2 + 9a(3siS 2 

We make some preliminary observations that enable us to compute the 
coefficients of the power series expansions of P(s\, s 2 , t) and Q(s±, s 2 ,t). Let 
(1*1,1x2) £ [0, l) 2 , and 9 £ [0, 1), and assume that 

ui + u 2 — 9u\u 2 < 1. (3.15) 

(It is clearly greater than zero.) Let 

1 00 

:= ]T (3.16) 



1 - Ui - U 2 + 6»MiM 2 j fc=0 

and 

00 

- log(l - Ul -u 2 + 9 Ul u 2 ) := £ C^wi^. (3.17) 

j,fc=0 

We give explicit expressions for Cj^ and Djj,. To begin we give several 
equalities that are easy to verify. 



Lemma 3.1 For u G [0, 1) 



1 =frrv 



(1 - u) q V « 
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p=0 



\P/ \Pj 



(3.19) 



Proof To get 

differentiate 1/(1 — u) q , n times, divide by n\ and set 
u = 0. □ 

We list the following equalities without proof. 
Lemma 3.2 

(p) ( P ) (p - l) 

VP/ fc w 

(k-l\ v (k\ 

\p-l) k\p) 

fk\ k-p+lf k \ 

\pj P \P - V 

Lemma 3.3 For < j < k 

D hk = j:{l-ey( 3 \( k \ (3.20) 



Also, C ,o = 0, C jfl = 1/j, C 0ik = l/k,j, k^O, and for 1 < j < k 
Proof Note that 

oo 

1 ^ = E( M i + M 2-^ lM2 ) n . (3.22) 

Writing (13.151) in the form u\ + w 2 — U\U2 + (1 — 6)uiU2 < 1 we see that it is 
equivalent to the statement that 

(1 ^ )MlU2 <1. (3.23) 



(l-« a )(l-« 2 ) 
We write 

1 



1 — U\ — U 2 + QU\U2 

13 



(3.24) 



(1-Ul)(l-U 2 ) 



(1 - 9)uiu 2 
(l- U i)(l-u 2 ), 



1 ™( (1 - 0H« 2 



1 - Wi)(l - M 2 ) \( X ~ _ M 2), 

;i - ey u \u p 2 



5 71 



£ (1 - u x )*+\l - u 2 )p+ 1 

Using this series we can determine Dj^. Since j < k it is clear that each 
of the first j + 1 terms in the series immediately above can contribute a term 
in u\u\. For a given < p < j we get (1 — 8) p times the coefficient of u\~ p in 
the power series expansion of 1/(1 — ui) p+1 and the coefficient of u 2 ~ p in the 
power series expansion of 1/(1 — u 2 ) p+1 . We see from Lemma [3. II that they 
are = (j) and = (J) respectively Thus we get (ET20j) . 

To obtain (I3.2ip we write 

— log(l — U\ — u 2 + 9uiu 2 ) (3.25) 
= - log((l - Ul )(l - u 2 ) - (1 - 6) Ul u 2 ) 

(1 - 9)uiu 2 \ 



log(l - Mi) - l0g(l - U 2 ) - log 1 



1 - «i)(l - u 2 ) 



h n h n ^ipV(l-«i)(l-«a); ' 

This gives us C^o and Co,a; and, similar to the computation of Dj t k, we can 
use the last series above to see that 

3(i-oyfj-i\{k-i\ 



P =i 



^ (l - /j - 1\ 

P+ 1 V P )\ p J 

p=0 



K p J \p J p + 1 jk 



This gives us (13.211) . □ 
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Set 

1 , 
* = ; 3.27 

By (13.121) we have that 

dt <t <dt + 2 } (3.28) 

for all t sufficiently large. 

Lemma 3.4 For all t sufficiently large, Pj,o(t) = oPjj, Po,k(t) — P k /k and 
for all 1 < j < k 

(See ([33]).) 

Proof Note that since < a, 0, 9 < 1, 

asi + (3s 2 — 6a(3siS2 < a + (3 — 6af3 (3.30) 

= a + p - a/3 + (1 - 0)a/3 

= -(l-«)(l-/3) + |^ + l + 0(l/t 3 ) 

ab 1 „ . . o, 

+ ^ + ! + 0(lA 3 ) 



dH 2 d 2 t 2 

Consequently, for all t sufficiently large 

< asi + /3s 2 - 6a/3sis 2 < 1. (3.31) 

Therefore, by (13. 17j) 

oo 

- log(l - cmi - /9m 2 + aP9u 1 u 2 ) = o? p k C jik u{u 2 . (3.32) 

j,fc=0 

The lemma now follows from Lemma 13.31 and (I3.27P . □ 

Set 

t = Vl + at + bt + dt 2 (3.33) 

and note that 

d l/2 t <t< d 1/2 t + 2 (3.34) 

for all t sufficiently large. 
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Lemma 3.5 For all t sufficiently large, and j, k > 1 

( 7 t 2 ( a -l)+t(2a-l)) 

Vi,oW = p " J (3-35) 

and 

g , i(() = . W-'HW-V , (3. 36) 
Furthermore, for all t sufficiently large and for all 1 < j < k 

Qj,k{t) (3.37) 

(--ft 2 f (1 - a)(l - 0) - (1 ~ P)P - {l - a)P + r^tll) 
V V J k J k J 



Proof It follows from (I3.3ip that for < Si, s 2 < 1 



=5 (asj + /3s 2 - ^a/5sis 2 ) n . (3.38) 



1 + aAi + b\ 2 + rfAiA 2 t 2 n (l 

Using this along with (13.141) we see that 

nl A (-ft 2 (l-s 1 )(l-s 2 )+t(2-s 1 -s 2 )) 

Q(s 1 ,s 2 ,t) = p (3.39) 

oo 

• 5Z(« s i + ^ s 2 - 6af3s 1 s 2 ) n 

n=0 

(-ft 2 + 2t- (-ft 2 + t) Sl - (-ft 2 + t)s 2 + -ft 2 Sl s 2 ) 



t 2 

oo 

Y,( as i + Ps2 - 0a/3s 1 s 2 ) n . 

n=0 



Consequently 



n m ((it 2 + 2t)a - (-ft 2 + 1)) 

Qj,o\t) = p a ' ( 3 - 40 ) 
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from which we get (13.35|) . and 

QoAt) = J2g±yg±ffl ^, (3.4!) 

from which we get (I3.36p . 

To obtain (I3.37P we use (I3.39p . and the terms of Dj jk defined in (13.161) . to 
see that 

Sm-J^T (3-42) 

= -7i 2 {D j>k a>P - Dj^a - Dj^ k f3 + ZVu-i) 
-t (2D jik af3 - D^ k _ x a - . 

Using (E32D, (B3QD and Lemma Owe see that 

Q jtk = ±t- 2 *( j )( k ) (3.43) 



p=0 



,2 / , _ k ~P oJ ~P , 3 -pk-p s 



-jt af3 — a — {3 



k j j k J 

-t ( 2a 0- a !^-pLjp\y 

Consequently, for all t sufficiently large, for all 1 < j < k 

-t Uap-a^-ptjlYj. 
Consider (13.441) . We write 

a k ~P J-P,j-P k -P ,o AK\ 

ap — a — : p 1 : — (3.45) 

k j J k 

n „vi m (1-qQp i P 2 

= (l-a)(l-0)-— — + jk> 
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2ap - a^-^- - (3 3 -y- = -a (l - /3 - |) - (3 ( 1 - a - ? j , (3.46) 



to obtain 

Q h k(t) ( 3 - 47 ) 
Note that for 1 < q < j 

Therefore, for each 1 < p < j 'we incorporate the term in p 2 /jk in (13.471) into 
the preceding term in the series in (13.471) to get (I3.37p . (Note that we can 
not add anything to the p = j term. The expression in (I3.37j) reflects this 
fact since = when p = j.) □ 



4 A sufficient condition for a vector in R 2 to 
be infinitely divisible. 

We present a sufficient condition for a random vector in R 2 to be to be 
infinitely divisible, and show how it simplifies the task of showing that ((G\ + 
c) 2 , (Cr2 + c) 2 ) is infinitely divisible. 

Lemma 4.1 Let if) : (-R+) 2 — > (-R+) 2 be a continuous function with if) (0,0) = 
1. Let s G [0, l] 2 and suppose that for allt > sufficiently large, \ogif)(t(l — 
si),t(l — sz)) has a power series expansion at s = given by 

oo 

<P(t;s u s 2 ) = W*K4 (4- 1 ) 

j,k=0 
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Suppose also that there exist an increasing sequence of finite subsets Mi C 
N 2 , i > 1, with U^i-M = N 2 , and a sequence t{ — > oo, i > 1, such that 
bj,k(ti) > for all (j, k) G M/(0, 0) and 



Km £ IM*i)l=0- 



(4.2) 



Then X 2 ) is the Laplace transform of an infinitely divisible random 

variable on (-R+) 2 . 

Proof It is clear from (14. 2 j) that the power series in (14. ip converges abso- 
lutely for all s G [0, l] 2 . Let 



i(U, si, s 2 ) = b ,o(ti) + b jtk (U)s{s2. 

j,feeJV</(o,o) 



(4.3) 



Set 



exp 



(ti;e" Al/ti ,e- A2/t 
We show that for each (Ai, A 2 ) G (-R+) 2 

lim (u;e- Xl/t \e- X2/t 



-Ai/ti —M/ti 



X 2 



(4.4) 
(4.5) 

As we point out in [9, page 565], (tf, e~ Xl ^ u ,e~ X2 ^ u ^ is the Laplace 
transform of a discrete measure. It then follows from the continuity theorem 
and the fact that ^(0,0) = 1 that ^(Ai,A2) is the Laplace transform of 
a random variable. Furthermore repeating this argument with (p^tf, Si, s 2 ) 
replaced by <f>i(ti] Si, s 2 )/n shows that -?/> 1//n (Ai, A 2 ) is the Laplace transform 
of a random variable. This shows that ip(X±, A2) is the Laplace transform of 
an infinitely divisible random variable on (R + ) 2 . 

Let 



^(1 - e- Al ^), U(l - e- x ^)) - ^(A l5 A 2 ) 



(4.6) 



Clearly Km™ 5* = 0. By (Oj) 



il>(ti(l-e- Xl / u ),ti(l 



)) - exp UiiU; e 



exp 6 ,o (U) + & i,fc(^) e 



(4.7) 



W.fc)^(o,o) 
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-exp l&o,o(**)+ & j,fc(^) e 
V 3,keNi/(p,o) 



-jAi/ti „-fcA 2 /t; 



where 



1 — exp 

i,^Mu(o,o) 



(l - exp ( - Y, hk(ti)e- jXl/ti e- kX2/U )) 
Note that by (jO]) 

lim €i = 0. 

i— >oo 

Therefore, by the triangle inequality, (14.61) and (14. 9 j) 



exp e 



-Ai/tj „-A 2 /ti- 



(4.8) 
(4.9) 
(4.10) 



< e^(ti(l - e- Al/il ), ti(l - e"^)) + 5*. 

Using ( 14. 6 ft we see that this is 

< e { (if;(Xi,X 2 ) + 6 i )+5 i . 

Thus we justify ( 14.51) and the paragraph following it. □ 

Remark 4.1 In P, Lemma 13.2.2] we present the well known result that the 
conclusion of Lemma [4.11 holds when \ogip(t(l — si),£(l — S2)) has a power 
series expansion at s = with all its coefficients, except for the coefficient of 
the constant term, are positive. Lemma I4TT1 is useful because it allows us to 
only verify this condition for a subset of these coefficients, (depending on t). 

The following lemma enables us to apply Lemma 14.11 

Lemma 4.2 For any C3 > there exists a constant B = -6(7, d, C3) for which 

tf t = {(j,k)\Jjk<Bt\ogt} (4.11) 



has the property that 



uniformly in |e| < C3. 



lim V \R jk (t,c)\ = 0, 

t — >oo — 



(4.12) 
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Remark 4.2 It follows from Lemmas 14.11 and 14.21 that in order to prove 
Theorem 11.41 we need that only show that we can find a Cq > 0, such that 

#j,fc0,c )>0 for all ^fjk<Btlogt, (4.13) 

for any constant B, for all t sufficiently large, (except for R Q (t)). 

Before proving Lemma 14.21 we establish the following bounds 

Lemma 4.3 

' ^ < ( ^ ] (4.14) 



\PJ V P 



and 



Proof It is clear that 



< J —. (4.16) 
pi pi 



Therefore to prove (14.31) we need only show that 



e 



v 



P ] -{-) > L ( 4 - 17 ) 

In [6J page 42 ], Feller shows that p\ {e/p) p is increasing in p. Since it is equal 
to e when p = 1, (and 1 when p = 0), we get (14.171) . 

The first inequality in (14.151) follows from (14.141) the next one is obtained 
by maximizing the middle term with respect to p. □ 

Proof of Lemma 14.21 By (13.81) and Lemmas 13.41 and 13.51 we see that for all 
t sufficiently large, for all 1 < j ' < k, 

\R j>k (t,c)\ < Ca>p k j2r 2p ^ Q (4.18) 

where C depends on c, 7 and d but not on j, k, t or t. Furthermore, C is 
bounded for all |c| < T, for any finite number T. (We also use the fact that 
lim^ooa/5 = 1.) 
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For any 8 > 0, for t sufficiently large, 

= at + dt 2 = _ l + bt < _ (f - J)6 < e -(i- a) 6/( dt ) 
f + at + bt + dt 2 l + at + bt + dt 2 ~ dt 

(4.19) 

and 

= frt + dt 2 = 1 _ l + at < 1 _ (1 - j)a < e _ ( i_j) a/(dt ) 
1 + at + 6t + c/t 2 1 + at + to + dt 2 ~ dt 

(4.20) 

Using these estimates along with (13. 28j) we see that for all t sufficiently 
large, for all 1 < j < k, 

\R jik (t,c)\ < c7 e -*M./(*) c -i(i-W(*) Q (4.21) 

uniformly in \c\ < c%. 

Suppose that y/Jk/(dt) = n. Then 

e -k(l-S)a/(dt) e -j(l-8)b/(dt) (4.22) 

= exp (-(1 - 8) (ajk/j + b^fjjPj n) 

< exp (-2(1 - 8)Vabn) 

' n 

where, for the inequality we take the minimum oia9+b/6 and for the equality 
we use the fact that ab — d+1. Combined with (I4.15P this shows that when 
\fjkj (dt) = n 



exp (-2(1 - 8) Vd + 1 



\Rjjk{t,c)\ < Cj exp (-2((1 - 8) Vd + T - l)n) . (4.23) 
Let A n = {(j, k)\n < </jk/(dt) < n + 1}. Then for any M 



E E E (4.24) 

jk/(dt)>M n=M (hk)&A n 



We note that the cardinality of A n is less than 

3nd 2 t 2 log((n + l)dt). (4.25) 
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This is because (j, k) G A n implies that 
n 2 d 2 t 2 (n + l) 2 d 2 t 2 

^-L < k < and j<(n+l)dt (4.26) 

J 3 

and summing on j we get (I4.25[) . 

It follows from (ETUI) and that 

£ |^ ifc (t,c)| (4.27) 

jk/{dt)>M 

oo 

< C(dt) 4 ^ n 4 exp (-2((l-5)Vd+T-l)n 

\4 



(dt) 4 



Clearly, there exists a constant 5 such that when M = B logi, this last term 
is o(l) as t — ► oo. □ 



5 Proof of Theorem 11.41 when (ci, C2) = (c, c) 
and EG1G2 > 

In this section we prove Theorem 11.41 in case 1. and for EG1G2 > 0, by 
establishing the positivity conditions on the coefficients Rj t k{t,c), (when 
EGiG 2 = 1), as discussed in Remark [4.21 We pass to the case EGiG 2 > 
on page UU 

To proceed we need several estimates of parameters we are dealing with 
as t — > 00. They follow from the definitions in fl3.10p - fl3.12p . 

Lemma 5.1 As t — > 00 

h 1 4- h 2 

1"« = 4-T^ + Or 3 ) (5.1) 



dt (dt) 
a 1 + a 



2 



l - p = it-JW + 0[t 



-3> 



v 1 "^ 1 -^) = ^ (^ +0 (* 



-4\ 
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a j = e -bj/(dt)+0(f/t*) 
f3 k = e -a.k/{dt)+0(k 2 /t 2 ) 

Also 

- (d + 2) 7 + d(a + b) = 2((d + 2) - (a + b)) = -2C (5.2) 
aj — d = (a — l) 2 
b-f-d = (&-1) 2 . 



Proof 



l + bt . , 

1 — a = ; — (5.3) 

1 + at + bt + dt 2 y ' 

l + bt 1 



dt 2 1 + a(dt)- 1 + b(dt)' 1 + d-H- 2 
= l±^(l- a (dt)- 1 -6(^)" 1 + 0(r 2 )) 

6 , d-b{a + b) t „ u _ 3 ^ 
~ dt + dH 2 + 1 > 

dt dH 2 +U[ ] 

The rest of the lemma follows similarly. □ 

Proof of Theorem 11.41 when c x = c 2 = c. To begin let note that it is 
easy to see from Lemma [3.41 that Pj,k(t) — for all < j, k < 00, with the 
exception of Po t o(t). This must be the case because exp(P(a, b, Ai, A2)) is the 
Laplace transform of an infinitely divisible random variable, as we remark 
following the proof of Lemma 12.11 
By $3M) 

n (A (lt 2 (a-l)+t(2a-l)) ._ x 

Qj,ow = p a ( 5 - 4 ) 



((- 7 6 + l)t + 7 (l + 6g)-26 + 0(l/t)) x 

= Or 

t 2 

((6- l) 2 t + 26- 7 (l + 6 2 ) + 0(l/t)) , x 



t 2 

0(l/t 2 ) I a J 



a 3 ' 



(it 
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Similarly 



Qo,k(t) = (^r^ + °( 1 /t 2 ))« fc - 1 . (5.5) 



Thus we see that there exists a t\ sufficiently large such that for all t > t\, 
Rj,o(t, c) and Ro,k(t, c) are both positive for all j, k > 1. 

We now examine Rj,k{t, c) for j A k > 1, j < k. We write 

j 

Rj,k(t, c) = Y^ RjApit, c )- (5-6) 

p=0 



Using (13371) and d3^29|) we see that 

Rj,k,p(t] c) 



(5.7) 



«/3t' 2p /A A\ 1 (j-p)(k-p) + cH^ (j\ fk\ 
t 2 \p)\p)p+l jk t 2 \p)\pj 

f- 7 t 2 ( (i - q)(i - /?) - (1 - ^ - (1 ~ Q)p + r 2j ~ vk -v 



j k j 

«-(--!)))■ 



+tla(l-/3-^+/3(l-«- P 



When p = we get 



ai-xp-x ~ 12 + p ^ 
(- 7 t 2 ((1 - a)(l - /?) + t~ 2 ) + t (a(l - /3) + 0(1 - a))) 

which is independent of j, fc. Using Lemma [5. II we see that 

--ft 2 ((1 - a)(l - 0) + r 2 ) + t (a(l -13) + (3{l - a)) (5.9) 

(d+2) a + b /l 
-7 + — r- + 



d? 1 d \t 
-(d + 2) 1 + d(a + b) | Q /l 

<i 2 
+ ' 



d 2 \t 
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Using this and Lemma 15.11 again we get 

Rj,k,o(t,c) 1 " 2c 2 (C/rf) + O (1/t) 



(5.10) 



aP- 1 P k - 1 d 2 t 2 
where the 0(1 /t) term is independent of j and k. 

We now simplify the expression of the other coefficients Rj y k, P (t,c), 1 < 
P < j- Set 



Rj,k,p(t, c) _ r-2p( J] (k\ ( 1 



c 



2 



a*" 1 /?*- 1 W VW U 2 * 2 

where 



FjApW = a $J^\ J k ( 5A2 ) 



and 

Aj,k, P (t) 



= - 7 t 2 f (1 - a)(l - ft) - (1 - {1 ~ a)p + t^tHtlP) 

V 3 k j k J 

+t(a(l-P-fy+P^-a-j\y (5.13) 
Using Lemma [5. II we have 
- 7 t 2 ( (1 - a)(l - 0) - (1 ~ ^ - {1 ~ a)p + r> j - pk - p ) 

V j j k J 

7 (d+1 ( bdt - i 1 + b<2 ))P (adt- (l + a 2 ))p j -pk-p\ 
d 2 V k J j k J 

^(-,, +2)+ M^ + H^>_g + o ( i) ( , 14) 
= -_ + 2 ) + b(dt - b)p + a(dt : a)p - 4) + o a 

d \ k j jkl \t 



and 



tlall-p-^+pll-a-Z)) (5.15) 
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'L-#-)-#-»)| +0 (i 

a \ j k J \t 

In ( 15.141) and (15.151) the expressions O (1/t) are not necessarily the same 
from line to line. Nevertheless, it is important to note that they are indepen- 
dent of p, j and k. That is there exists an M > such that all terms given 
as O (1/t) in (l5~T4"D and (l5~T5|) satisfy 

_ ^ < o fl) < E. (5.16) 



This is easy to see since the O (1/t) terms, in addition to depending on t, 
depend on a, 6, p/j and p/j < p/k < 1. 
Using (jSTHj) . (g35D we have 

W) = i(-( d + 2 ) + bk ^ + a ^^-^) < 5 - 17 > 

+ lfa + t - P ' <i '.-°»-^- t )Uo^ 1 



<i I j k J \t 

-(d + 2) 1 + d(a + b) 
d 2 

(7a - d)p(dt - a) (76 - d)p(dt - 6) _ 7^ /l 
id 2 fed 2 jkd 2 \t 

-2( + (a - l) 2 p(rft - a) + (b-l) 2 p{dt-b) 7P 2 | Q /I 



d 2 jd 2 fed 2 jW 2 Vt 

where, for the final equality we use (15.21) . 
Note that 

(a-l) 2 p(dt-a) (6 - l) 2 p(rft - b) 
B jtKp [t) := + (5.18) 

> 2p[(a- 1)(6- 1)|- 

fjk 



<J{dt- 


a)(dt — b) 






- a)(dt — b) 
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(For the inequality use a 2 +(3 2 > 2a (5.) Therefore since ( = a+b—(d+2) > 0, 



A ^ (t) * T 2 [ p ^ 1 J c -^ + U 

- 2^ + 0(1)) v v +0 m (5 _ 19) 



d? \ ^/Jk J d 2 jk \t 

2 (** ( 1+0 ( 1 -)-^ T )- 1 )( + 0^ 



d 2 \VJk V V*/ (VJkdt J J \t 

Thus we see that there exists a function e<, depending only on a and 6 
such that 

^ (i >4(S (1 - e,, - 1 ) c+0 ©' (5 - 20) 

where 

lim e t = 0, (5.21) 

t^oo 

and, as we point out above the O (1/t) is independent of p, j and k. 

Remark 5.1 We interrupt this proof to make some comments which may 
be helpful in understanding what is going on. Note that if 

V jk < 1 - e for some e > (5.22) 



dt 
then 

1 _ 2r 2 (C /d) 

Rj,k{t, c) > R jA0 {t, c) > (1 - 6) d2 £' } - a?- x p h - x as t - oo (5.23) 

for all 5 > 0. This follows from ( 15.101) and (15.201) since when (15. 22ft holds 

^( t )>^(^|-l)c + o(i)>0, (5.24) 

for all p > 1, for all t is sufficiently large. Consequently when ( 15.221) holds 
Rj,k(t, c) > for all t is sufficiently large when 

c 2 < A (5.25) 

(Here we also use (I5.2ip .) 

(When C < 0, (15T231) shows that Rj, k (t, c) > for all c G i? 1 . This is what 
we expect. (See the paragraph containing ( 12. 12ft .) 
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We use the next two lemmas to complete the proof of Theorem 11.41 in 
case 1. 

Lemma 5.2 For any N £ R + , we can find Co > and t co < oo such that 
for all t > t Co 

R jAp (t,c)>0 (5.26) 
for all |c| < Co and all p, j and k for which jk/t < N Q . 

Proof This follows from (I5.1(jp when p = 0. Therefore, we can take p > 1. 

We first show that for any N G R + , Rj^, P {t) > when ^/Jk = Ndt, for 
all t sufficiently large. By Remark 15.11 we can assume that N > 1 — e. It 
follows from flo^Ol) that 

^W*) ^ 4 (tpO- -et)-l)( + (]) , (5.27) 



d? \N y ' J" \t, 
where e t satisfies (I5.2ip . Therefore when p > AN for any A > 1, A, jfc p (t) > 0, 
and hence Rj^,p(t, c) > 0, for all t sufficiently large. 
Now suppose that 

p < AN. (5.28) 

Since ^fjk = Ndt we see that 

Wk^ = ° (1/t2) ' (5 ' 29) 

where the 0(l/t) term is independent of p, j and k. 
Note that by <\5A$ and fl5TT2|) 

k^^ ( Ww +0 ^ (5 - 30) 

Therefore, if in addition to fl 5 . 2 8 1) we also have AiV < j/2, so that p < j/2, 
we see by (joTTT!) . fl5TT7j) and (15T291 that 

^^(*) + ^ A j,k, P (t) ( 5 - 31 ) 
> J_ / (j-p)(fc-p) _ c^2C\ Q / 1_ 



dH 2 V (p + l)jfc J \t 3 

- d 2 t 2 v4(p + i) d / u 3 

1 f_L c^2C\ / 1_ 
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Therefore we can obtain Rj^ p {t) > by taking 

c 2 < (5.32) 
~ 16AWC 

for some A' > A. 

Now suppose that AA" > j/2. In this case we use (15.171) to see that 

2C (a — l) 2 pt , . , , 
A jik>p (t) >~^+ 1 J. P + 0(l/t). (5.33) 

It is easy to see that the right-hand side of (15.331) is greater than zero for all 
t sufficiently large since 

{a ~ 1)2pt >Tl^(*-^- (5-34) 



2dj ~ 4dNA 

Thus we see that for any fixed N, Rj^, p {t) > for all t sufficiently large. 

Since the 0(l/t) terms are independent of p, j and k this analysis works 
for all j and k satisfying (15.261) . and all 1 < p < j as long as (15.321) holds 
with A^ replaced by N . □ 

Lemma 5.3 For all N and B G R + we can find a c' > and t c > < oo such 
that for all t > t c i 

R j , k , p {t,c)>0 (5.35) 
for all \c\ < c' and all < p < j < k for which 



N t < ^jk < Btlogt. (5.36) 

(The value of N in Lemmas 15.21 and 15.31 can be taken as we wish. It will 
be assigned in the proof of this lemma.) 

Proof By adjusting A" and B we can replace (15.361) by the condition 



N t < ^jk < Btlogt. (5.37) 
Using (15.171) , we see that if j < pt 

-2C (o - lfp{dt - a) (b - Ifpjdt - b) 
]Ap[) " rf 2 + jd? + kd* 
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1T>2 • O (-) (5.38) 



jkd 2 \t 

-2C (o - Ifpjdt - a) 7 /l 

- rf 2 jd 2 rf 2 U 

-2C (a — l) 2 7 ^/l 
> — - + - + O I - 

~ d 2 2pd d 2 \t 

Clearly, there exists a p > 0, independent of j and k such that this term is 
positive. Thus we can assume that 

j > pl (5.39) 

Furthermore, when \fjkjt = N, it follows from (13. 28f) that we can write 

( EZEgp as 

W> 4 (& (l + O (j) - ^) - l) C + O (i) . (5.40) 

Let 5tv = (10 log iV/iV) 1/2 • Clearly, if p > (1 + 5jv)-W, the right-hand side of 
(I5.40p is positive for all t sufficiently large. Therefore, when y/JTt/t = N, we 
may assume that 

P<(1 + 5 N )N. (5.41) 

(The value chosen for 5^ simplifies calculations made later in this proof.) 
In addition we can also assume that 

P > Po (5-42) 

for any finite po, since if p < po 

F jAp (t) > F j,k, P0 (t) > c 2 A jAp (t). (5.43) 
for all c > sufficiently small. 

We use the next lemma in the proof of Lemma 15.31 

Lemma 5.4 For j < k, with p and j large andp/j small 



\PJ \Pj 2vrp \ p 2 



(5.44) 



exp (-iL(p - 1) - ^(p - 1) + 0(p 3 /j 2 )) (l + 0(p- 1 ] 
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When t G R + is large and y/Jk/t = N, under assumptions \5. 39\) and \5.J^1 ) 

M3(:k(t)Vo( P -'>). (-> 

Proof By Stirlings's formula for integers q, 

q\ = V2^q q+1/2 e- q (l + 0(g -1 )) . (5.46) 
Therefore, since j is large and p/j is small, terms of the form 

(i + or 1 )) 

(l + 0( J5 - 1 ))(l + 0((j-p)- 1 )) 
Using this we see that 

(A j- 



l + 0{p- 1 )). (5.47) 



(5.48) 



1 ,-i+i/a 

f i + o(p- v 



2 n (j _ p)(j-P+l/2) 



V 72 ^ V p / (j - p)(l'-P+V2) 



p 



V^pypJ (1-2 

V 



(j-p+1/2) 



l + 0(p' 



-1 



\/2np \p J 
1 ftV 



1 / e j ' 



yfhvp \ p 



e -(i-p+i/2)io g (i-p/i) A + 
e 0-p+i/2)(p/j+ P 2 /(2i 2 )+o(p 3 /i 3 )) _|_ 0(p -1 )) 

9 

e (-p 2 /(2j)+p/(2j)+o( P 3 /i 2 )) + 0(p -1 )) 



e (-p(p-i)/(2i)+o(p 3 /j 2 )) A + Oip- 1 )) . 



yfhxp \ p 

Since this also holds with j replaced by k we get (I5.44p . 

To get (15.451) we multiply each side of (15.441) by t~ 2p and substitute for 
y/Jk/t = dN and use the fact that under the assumptions ( 15.391) and ( 15.411) . 

£ < t. < ( 1 + ^) 2Ar2 



- 2 <-< J • (5-49) 
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Consequently, for all t sufficiently large 

exp (-^(p - 1) - A(p _ i) + 0( P 3 /J 2 )) =1 + (N 2 /t) . (5.50) 

□ 

Proof of Lemma 15.31 continued We show that under the assumptions 
(I539|) and (ETlIj) . when vW* = for N < N < Blogt, for any < B < 
oo, and t is sufficiently large, 

V X T / . \ OAT -L 



T E W*)>^ M (5.51) 



ai- 1 ^- 1 p f^ iV3/2 
for some C > 0, independent of iV, and 



a ,-i^-r E Wt)>-^p (5-52) 

for some D < oo, independent of AT. If (I5.5ip and (15.521) hold, we can find a 
c > such that for all c 2 < Cq 

(1+5jv)JV 

E i2i,fc jP (*,c)>0 (5.53) 

P=P0 

for all i sufficiently large. Since we have already established that Aj^, P {t) > 0, 
when p < Po and p > (1 + Sn)N, this completes the proof of Lemma I5T31 

Thus it only remains to prove (15.511) and (15.521) . We do (I5.5ip first. It is 
considerably easier than (15.521) . By Lemma [3.41 and (15.451) 

p (1+S N )N 



T E P JA P (t) (5-54) 

P=P0 

l fj\f k \(j-p)(k-p) 



P—Po 



t 2 P\pJ\pJ (p + l)jk 



P t; t^\p)\p)p 

>C (1+ ^ )N }_ (e_N\ 2p 
P = P0 P 2 \ P J 
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In order to calculate this last sum we consider the function 

2y 



]_(eNy V J- 

y rn [ 

for m > and y > 2. We have 



f >"M = ^[Y) =^** N - V * V) (5-55) 



-m 



f m {y) = ^_- + 2(l + log7V-logy)-2j/(y) (5.56) 
+ 2(log7V-logy) /(j/). 



v y 

This has a unique root y m where 



TJl 

\ogy m + — = logN. (5.57) 
2y m 



(Clearly, y — N). Let y m = JV(1 + e m ). Then 



/77 

log(1+f "' ) + 2iv(rT^= - < 5 - 58 ' 



Consequently 
which implies that 



777 

e m = -^ + 0(N- 2 ), (5.59) 

Tfl 

y m = N-- + 0(1/N). (5.60) 

Making use of the fact that (—m/y m + 2(1 + log iV — \ogy m ) — 2) = 0, we 
see that 

( m 2 \ 

/m(!/m) = — - — )f(ym)<0. (5.61) 
\ Urn Dm J 

Therefore 

sup f m (y) = f m (y m ) < — ^e 2Af . (5.62) 

y>2 (A/ — m) m 

We also note that since —m/y — 2 logy is increasing for y > m/2, f' m (y) is 

positive for m/2 < y < y m and negative for y > y m . Consequently, f m (y) is 
unimodal. 
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Now consider the last line of (I5.54p . The function being summed is i^G )- 
The above discussion shows that this function is unimodal, with a maximum 
at, at most, two points at which it is less than 2e 2N /N 2 . Consequently, to 
obtain (15.511) we can replace the sum in the last line of (15.541) by an integral 
and show that 

Making the change of variables r = xN we have 

1 rl+&N l / e \ 2xN 

N Jp /N X 2 \x) 

Recall that N < N < 2 logt, and that we can take N as large as we want, 
(but fixed and independent of t), and that 5n = (10 log N/N) 1 / 2 . Therefore 



1 /-l+ao log TV/TV) 1 / 2 1 / e 



2xN 



h > - / - - dx (5.64) 

N Jl-(101ogAT/AT) 1 /2 x 2 \xj 
I |.l+(101ogAf/AT) 1 /2 , e \2xN 

~ 2N 7l-(101ogAr/Af) 1 /2 \xj 



We write 



e \ 2xN 



exp(2xN{l -logx)). (5.65) 



Set x = 1 + y and note that for \y\ small 

x(l-logar) = (l+y)(l-log(l + ?/)) (5.66) 

= (1+ 4 1 -| ( - irI ?) 
= 1+s _f ( _ ir .(l±«0!r 

n=l 71 
°° f 1 I 



n=2 

oo n 

i-EHr-r 



n n — 1 



When 1 2/ 1 < (10 log N/N) 1 ' 2 , so that \y\ 3 N « 1, this shows that 



2xN 



e 2N e -y-N (l + (\y\ S N)). 
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It follows from this that when we make the change of variables x = 1 + y in 
( IBTMj) we get 

e 2N ,(101ogJV/iV)l/2 

h> — e~ yN dy (5.68) 

_ 2iV 7-(ioio g iv/7V)i/2 y v ; 

e ™ ,(20 log V) 1 /^ 



> : / e -«r/a ^ 

2^^3/2 i-(201ogiV)i/2 



Since 



/ e~ M /2 du < N- 10 , (5.69) 

J (20 log TV) 1 / 2 

we see that (15.631) follows. Thus we have established ( 15.511) . 
Before proceeding to the proof of (15.521) we note that 

(1+8 N )N , / A /,\ 2N 

To prove this we use (I5.45P and the same argument that enables us to move 
from a sum to an integral that is given in (I5.54l) - (l5.63p . except that we use 
(I5.62p with m = l. We continue and then use (15.671) to get 



P=P0 



- Jpo u\u) \ N 



e 2N 
2N 1 / 2 ' 



We now obtain (15321) . When y/Jk = tN, by (I51T1) and fl5^0|) 

r- 



T Q j>k ,p(t) (5-72) 



t 2p \v \p 



2 P\( k \h(L 



i)c + (J 
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By f)5.7ip we see that 
Therefore, to obtain (I5.52p . it suffices to show that for some D < oo 

Here we use the fact that N < B logi for some < B < oo. 

Remark 5.2 Since the proof of (I5.74p is rather delicate we make some 
heuristic comments to explain how we proceed. When yfjk = tN the term 

^ _2P (p)(p)' as a ^ unc ^^ on °f Pi i s exp(2iV) times values that are sort of nor- 
mally distributed with mean p = N, and, roughly speaking, 

'T i(!)Q~Ce»^, (5.75) 



P=PO 



t 2p \P \P N 1 / 2 ' 



for all t sufficiently large. (In fact the upper bound is given in ( 15.711) .) This 
is too large to enable us to get (I5.74p so we must make use of the factors 
(jj — 1J, which is an odd function with respect to p = N, to get the can- 
cellations that allow us to obtain (I5.74p . However, because we are canceling 
terms, we must take account of the error in Stirling's approximation; (see 
(I5.46P ). To do this we need to show that the estimate in (15.741) remains the 
same even when we eliminate the terms in the summand that are not close 
to N. 

Proof of Lemma 15.31 continued Note that by (I5.45P 

jV(l-jV-V4) i NQ.-N-V*) 

p=p t 2p \Pj \P/ p=p 

The fact that f m (y) is unimodal on y > m/2 implies that fi(p) is increasing 
on the interval \po,N(l — N~ 1 ^)]. Therefore 

^ _ (™) < CNMNil-N-^)) (5.77) 
p=po P\P J 
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\ 2Af(l-Af" 1 /4) 

< c 



( ^ e 2Ar(l-Ar- 1 /4)(i_i g(i_Ar-i/4)) 
< ^y e 2AT(l-AT- 1 / 4 )(l+Af- 1 /4) _ q 2N-2N 1 / 3 

Let 5' N = N^ 1 / 4 . The argument immediately above shows that to prove 
(I5.74p . it suffices to show that 

By ([525]) 

-OI.Kt(H)H(3) <"* 

Using ( 15.7ip together with the fact that since p > (1 — S' N )N, 1/p < 1/N, 
we see that 

(1+S N )N / e pj\ 2 P{l l\ fl\ (~L+Sn)N / eN \ 2 P i /]\ 



1 



Therefore, to obtain (I5.78P that it suffices to show that 

E U"? - _Ce ^ (5 ' 80) 



p=(l-5' )N 



In a minor modification of the analysis of f m (y), we write 

%) == ( eJ f) V {^- l i)= eXP (% (1 + l ° gN - l0gV)) {if ~ I) ■ 

(5.81) 

Therefore 

h'(y) = ((2 (1 + log iV - log „) - 2) (1 - i) + i) ^ * . (5.82) 
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Let y = (1 + u)N. Then h'(y) = when 



2u;log(l + ^) + 1 =0. (5.83) 
N(l + to) 



This equation is satisfied when 



<2N \N 
Note that when y — (1 + u)N 

(eN\ y / e \ (i+"0* 



w = ±-^= + ( 4 ) . (5.84) 



< e N , (5.85) 



V y J vi + w> 

because (e/x) x is maximized when x — 1. Therefore 

2 2/ / 1 i \ J2N 



'^Vf-- i )<— — (5-86) 



from which we get 



2iV 



sup |%)|<C . (5.87) 

i<y<(i+s N )N \N^ 2 J 

It is easy to see that h(y) is negative for 1 < y < N and that it decreases 
to its minimum value at N(l — u) and then increases to zero at y — N. 
It then increases to its maximum value at N(l + u) and then decreases for 
N(l + oj) < y < (1 + 8n)N. Consequently the difference between 

(1+5jv)JV r (l+8 N )N 

Y h(p) and / Up) dp (5.88) 

differs by at most 4max 1 < p <( 1+ j JV ) JV |%)|. Since this is 0{e 2N /N 3 / 2 ) by 
f )5.87p . and we are only trying to obtain (15.801) . we can neglect this dis- 
crepancy. Therefore to obtain fl 5 . 8 1) we need only show that 

(i+s N )N i (eN\ 2p fp \ , . e 2N , 
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Under the change of variables p = xN the integral in (15.891) is equal to 

-l+<5iv 1 / e \ 2xN 

(x-1) dx. (5.90) 



1-S' N X \x 



By (I5.66P , in which x — 1 + y , 



c \ 2xN 

£) = e 2N ^N+yZN/Z+Oiy^N (5.91) 
XJ 

= e 2N e-y 2N (l + ^ + 0(y 4 )iv) . 

Therefore, with the change of variables x = 1 + y we write the integral in 
flCTD as 

e 2 " -^e- 2 - f 1 + tE. + tf) N ) dy. (5.92) 



l-s' N i + y V 3 

We use (1 + y)~ l = (1 — y + y 2 — y 3 + 0(y 4 )) to write 



l + y V 3 



1 + V- + 0(y 4 )iV (5.93) 



?/ 4 7V 

y-y 2 + y 3 + V — - y' + Otf)N. 



Using this we see that (15.921) 

= e 2N f_ N &/ (y-y 2 + y* + y -^f-y A + 0(y 5 )iv) dy. (5.94) 
Recall that 8 N = (10 hgN/N) 1 / 2 and 5' N = iV -1 / 4 . Since 



e 2N 



e- y ^d?/ < -— (5.95) 

(101ogJV/JV)V2 N 10 



and 

„_Ar-i/4 



e^ v / e"^ v dy < e 2 "~ N ' , (5.96) 



errors we can ignore in obtaining (15.521) . we can simplify matters by replacing 
the integral in (15.941) by 

e 2N jT e~y 2N (y-y 2 + y*+ V -^--y 4 + 0(y 5 )iv) dy (5.97) 
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-e 2N / " e-y' N (y 2 - y -^ + y 4 + 0(j?)n) dy 



00 2, 



00 
00 



3 

-e 2N I " e~ yI ( — V — + + 0(y°)N- 2 \ dy 




iV 2 " 



Since 





(5.98) 



00 

1 

= 4' 

we obtain (15.89H . □ 

Proof of Theorem 11.41 when Ci = c 2 = c and EG1G2 > concluded 

Consider the Gaussian random variable (G1/7, G2/7) where 7 = (PG1G2) 1 / 2 . 
This random variable has covariance V in (12. ip . By Lemma [5.31 there exists 
a Cq > such that (G1/7 + c, G2/7 + c) has infinitely divisible squares for 
all |c| < Cq. Let c be the supremum of the c for which this holds. Since, by 
hypothesis, ( 11. 12ft does not hold, c is finite. Therefore, (G1/7 + c, G2/7 + c) 
has infinitely divisible squares for all |c| < c and not for any c or which 
|c| > c. Translating this into the notation used in Theorem 11.41 we have 
(G1/7 + ca, G2/7 + ecu) has infinitely divisible squares for all |a| < c/c and 
not for any | a \ for which | a \ > c/c. 

Therefore, to complete the proof of Theorem 11.41 when c\ = C2 = c and 
EG1G2 > we need only show that (Gi/j+c, G2/7+C) has infinitely divisible 
squares for |c| = c. Consider the Laplace transform of (G1/7 + c, G2/7 + c) 
in ( 13. ip . Since it only depends on c 2 we can simplify the notation by taking 
c > 0. Let c m I c. Abbreviate the third line of (13.11) by exp (P + c 2 Q). Thus 
exp (P + c 2 m Q) is the Laplace transform of an infinitely divisible random 
variable. Therefore, for each t > the power series expansion of P + c 2 m Q in 
Si and S2 has positive coefficients, except for the constant term. Thus if we 
write 

j,k j,k 
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we see that a^k + c^bj^ > for each (j,k) ^ (0,0). Letting c m j c we 
therefore have dj^ + c^bj^ > for each (j, k) ^ (0,0). This shows that 
exp ((P + <?Q)) is the Laplace transform of an infinitely divisible random 
variable. □ 

Remark 5.3 In the remainder of this paper we continue to prove Theorem 
ll.4l for all Ci, c 2 and arbitrary covariance EG\G 2 . In each immediately 
above, because fll.l2f) does not hold, there exists a d < oo such that {G\ + 
cc%,G 2 + cc 2 ) does not have infinitely divisible squares for all c such that 
|c| > d. Therefore, if we can show that there exists some c ^ for which 
both 

(Gi + cci, G 2 + cc 2 ) and (G\ — cci, G 2 — cc 2 ) (5.99) 

have infinitely divisible squares, we can use the arguments in the preceding 
three paragraphs to show that there exists a critical point c such that {G\ + 
cci, G 2 + cc 2 ) has infinitely divisible squares for all |c| < c and not for |c| > c. 
Consequently, in the remainder of this paper, in which we consider different 
cases of c\, c 2 and arbitrary covariance EG\G 2 we will only show that (15.991) 
holds for some c ^ 0. 



6 Proof of Theorem 11.41 when (ci, C2) = (c, ±c) 



We first assume that EG\G 2 > and that (ci, c 2 ) = (c, — c). In this case we 
have 

1 ( 2 / pA!A 2 + Ai + A 2 \\ 

" (H(a, b, A l5 A 2 ))V2 exp ^ \ 2H (a, 6, A a , A 2 ) ) ) ' 

where p = a + b + 2. This is exactly the same as (13.11) except that 7 is 
replaced by p. We now trace the proof in Sections [3H5] and see what changes. 
Obviously much remains the same. In particular the power series P is un- 
changed. The basic expression for Q in f!3.14p is essentially the same except 
that 7 is replaced by p. Thus Lemma 13.51 is also essentially the same except 
that 7 is replaced by p. 

The analysis in Section @] only uses the fact that 7 < 00, and since p < 00, 
Lemma [4. II also holds in this case. 
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In going through Section Owe see the coefficients of Q change, but they 
still lead to essentially the same inequalities that allow us to complete the 
proof. In place of (I5.2p we have 

-3p + a + b = -2(3 + (a + 6)) := -2( (6.2) 
ap-1 = (a + 1) 2 
bp -I = (b + 1) 2 . 

Using this in (15. 4p and (15. 5p . with 7 replaced by p, we get 

Q jfi (t) = (^^ + 0(lA 2 ))^-\ (6.3) 



and 



Qo, k (t) = ( i -^r- + 0{1/t2) ) ak ~ 1 - (6 - 4) 



We also see that we get (15.71) with 7 replaced by p and consequently, in 
place of (15.101) . we get 

Rj,k fi {t,c) = l-2c 2 ((/d) + Q(t- 1 ) 

Of course the key term in the proof is the analogue of Aj t f. tP (t). We get 
the third line of (I5.17P with 7 replaced by p, which by (I6.2p leads to (15.191) 
with ( replaced by ( and 7 replaced by p. Therefore, all the subsequent 
lower bounds for Aj^#{t) that are in Section [5] hold when ( is replaced by 
(. In the proof of (15.521) in Section \5\ the only property of ( that is used is 
that is is positive. Since ( is also positive the same argument completes the 
proof of Lemma T5.3I and consequently, by Remark 15.31 of Theorem 11.41 when 
EGiG 2 > and (ci, c 2 ) = (c, -c). 

When EGiG 2 < and (ci, c 2 ) = (c, — c) we note that 

((d + c) 2 , (G 2 - c) 2 ) ^ ((d + c) 2 , (-G 2 + c) 2 ). (6.6) 

Now EG\(—G2) > and we are in the case proved on page [4X1 Therefore, 
by Remark I5.3[ Theorem 11.41 holds in this case. 

Finally when EGiG 2 < and (ci, c 2 ) = (c, c) we note that 

((d + c) 2 , (G 2 + c) 2 ) ^ ((d + c) 2 , (-G 2 - c) 2 ). (6.7) 

Now EG\{— G 2 ) > and we are in the case proved in the beginning of this 
section. □ 
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7 Proof of Theorem 11.41 when (ci, c 2 ) = (c, 0) 



We first assume that EG1G2 > 0. In this case we have 

E G ( e ~ (Al(Gl+c)2+AaG 2 )/2 ) ( 7 _-Q 

1 / J 5A 1 A 2 + A 1 \\ 

(if(a,6,A 1 ,A 2 )) 1 /2 exp ^ C V x 2^(a,6,A 1 ,A 2 ) y / y / ' 

The term in the numerator of the exponential lacks the A2 that is present 
in (13.1 1) and (16. ip . Therefore, the formulas for the coefficients of the power 
series for the analogue of Q, which we denote by Q, are different. It is easy 
to see that in place of (13.391) we get 

g (si , S2 , t ) = _ (*'(! -*)(! -*)+«(!-..)) (7 . 2) 

00 

• 5I( as i + 0*2 - 6>a/3sis 2 ) n 

n=0 

(6t 2 + t - (bt 2 + *)s x - 6t 2 s 2 + 6t 2 sis 2 ) 
P 

00 

• + - 0a/3s]S 2 ) n . 

n=0 

Using this, in place of Lemma 13.51 we get 
Lemma 7.1 For allt sufficiently large, and j,k > 1 

~ (tt 2 +*)(!- a) 

g ii0 (t) = p a J (7.3) 

and 

g^W = tfd-^+V - (7.4) 
Furthermore, for all t sufficiently large and for all 1 < j < k 

Qj,k(t) (7.5) 



t 2 



"0 vv VP/ 



-6t 2 f (1 - a)(l - /?) - (1 ~ P)P - (1 " Q)P ' - 2 j ~ P * " P> 



j fc J fc 

, ,))■ 



( 1 - a '' 
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The analysis in Section @] only uses the fact that 7 < 00. Since b < 00, 
Lemma [4. II also holds in this case. 

In going through Section [5] we see the coefficients of Q change, but they 
still lead to similar inequalities that allow us to complete the proof. Using 
(H3D, (J73D and flEQ) we get 

Quo® = [§i + 0(1/*)) ^ ( 7 - 6 ) 

and 

Qo, k (t)= (^T + °( 1 /*))/5 fc - 1 , (7-7) 

since ab — d + 1. 

We next consider he analogue of ( 15. 61) which we denote by Rj^it). We 
see that in computing this the first two lines of the analogue of ( 15.71) remain 
unchanged, except for replacing c by 1. The last two lines of ( 15.71) are now 

{-be ( (i - tt )(i - ft - {1 ~ P)P - {1 - a)p + r^tll\ 

\ \ J k j k J 

+^ll-a-?jj. (7.8) 

Therefore, in place of (15.101) . we get 

R lA0 (t) _ l-2c 2 (b/d) + 0(t- 1 ) 



Using (I5.14p , with 7 replaced by b and Lemma 15.11 we see that (j7.8p 



(7.9) 



4 ( (ab + 1) - hv{d \ ~ h) - ap{dt ~ a) + 41 (7-10) 
d \ k j jk J 

2b pjdt - a) b 2 p(dt - b) _bp^ ( 1 \ _ 
d 2 d 2 j d 2 k jk \tj 



-d 2 \ p ^jk d 2 jk \t 
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Comparing this inequality to the first line of (15.191) we see that we have 
exactly what we need to complete the proof in this case. The rest of the 
argument in Section [5] only uses the fact that £ > 0. It is now replaced by 
b > 0. Thus we get Lemma 15.31 and, by Remark 15.31 Theorem 11.41 when 
(ci,C2) = (c, 0) and EG1G2 > 0. However this proof holds for c positive or 
negative, so if EG1G2 < 0, we simply note that 

(Gl(G 2 + c) 2 ) l ^(Gl,(-G 2 -c) 2 ). (7.12) 

Since EGi(— G 2 ) > we are in the case just proved so, by Remark 15. 3\ 
Theorem 11.41 holds in this case also. □ 



8 Proof of Theorem 11.41 for (ci, c 2 ) <E R 1 x R 1 



It is simple to complete the proof from the results already obtained. Suppose 
neither c\ nor c 2 are equal to zero. Then, clearly, 

((G 1 + cc 1 ) 2 ,(G' 2 + cc 2 ) 2 ) (8.1) 

is infinitely divisible, if and only if 

((G 1 /c 1 + c) 2 ,(G 2 /c 2 + c) 2 ) (8.2) 

is infinitely divisible. We have already shown that there exists a critical point 
c > such that 

{{Gi/ Ci + c) 2 , (G 2 / c 2 + c) 2 ) (8.3) 

is infinitely divisible for all |c| < c and not for |c| > c. Consequently c is also 
a critical point for the infinite divisibility of 

((G 1 + c 1 c) 2 ,(G 2 + c 2 c) 2 ). (8.4) 

If Ci = we repeat this argument for 

(Gl,(G 2 + cc 2 ) 2 ). (8.5) 

□ 



46 



References 

1. Bapat, R. (1989). Infinite divisibility of multivariate gamma distributions 
and M-matrices. Sankhya, 51, 73-78. 

2. Eisenbaum, N. (2003). On the infinite divisibility of squared Gaussian 
processes. Prob. Theory Related Fields, 125, 381-392. 

3. Eisenbaum, N. (2005). A connection between Gaussian processes and 
Markov processes. Elec. J. Probab., 10, 202-215. 

4. Eisenbaum, N. and Kaspi, H. (2006). A characterization of the infinitely 
divisible squared Gaussian process. Ann. Probab., 34- 

5. Eisenbaum, N., Kaspi, H., Marcus, M. B., Rosen, J., and Shi, Z. (2000). 
A Ray-Knight theorem for symmetric Markov processes. Ann. Probab., 
28, 1781-1796. 

6. W. Feller (1950). An Introduction to Probability Theory and its Applica- 
tions, Vol. I, John Wiley, New York. 

7. Griffiths, R. C. (1984). Characterizations of infinitely divisible multivari- 
ate gamma distributions. Jour. Multivar. Anal, 15, 12-20. 

8. M. B. Marcus and J. Rosen (2006). Infinite divisibility of Gaussian 
Squares with non-zero means, submitted for publication. 

9. M. B. Marcus and J. Rosen (2006). Markov Processes, Gaussian Pro- 
cesses and Local Times, Cambridge University Press, New York. 



47 



